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ABSTRACT 

Context. The solar system was most likely bom in a star cluster containing at least 1000 stars. It is highly probable that this cluster 
environment influenced various properties of the solar system like its chemical composition, size and the orbital parameters of some 
of its constituting bodies. 

Aims. In the Milky Way, clusters with more than 2000 stars only form in two types - starburst clusters and leaky clusters - each 
following a unique temporal development in the mass-radius plane. The aim is here to determine the encounter probability in the 
range relevant to solar system formation for starburst or leaky cluster environments as a function of cluster age. 
Methods. N-body methods are used to investigate the cluster dynamics and the effect of gravitational interactions between cluster 
members on young solar-type stars surrounded by discs. 

Results. Using the now available knowledge of the cluster density at a given cluster age it is demonstrated that in starburst clusters 
the central densities over the first 5Myr are so high (initially > 10^ MqPC"^) that hardly any discs with solar system building potential 
would survive this phase. This makes a starburst clusters an unlikely environment for the formation of our solar system. Instead it 
is highly probable that the solar system formed in a leaky cluster (often classified as OB association). It is demonstrated that an 
encounter determining the characteristic properties existing in our solar systems most likely happened very early on (< 2Myr) in its 
formation history and that after 5Myr the likelihood of a solar-type star experiencing such an encounter in a leaky cluster is negligible 
even if it was still part of the bound remnant. This explains why the solar system could develop and maintain its high circularity later 
in its development. 
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1. Introduction 

The average stellar density in the solar neighbourhood is 0.112 
stars pc"-', meaning that the Sun resides in relative isolation. 
This might not have been the case at the time of its formation. 
Most stars form within groups or clusters of stars (Lada & Lada 
2003, Pon-as et al. 2003, Evans et al. 2009). The majority of 
these clusters and associations dissolve within some tens of Myr 
(Lada 2006, Portegies-Zwart et al. 2010). Even the long-lived 
open clusters do not exceed a lifetime of some hundreds of Myr 
At the current epoch the Solar System is about 4.6 Gyr old. This 
means if the Sun formed in a cluster environment, the solar birth 
cluster would have dissipated a long time ago (Portegies-Zwart 
2009). 

Several characteristics of the today's solar system give clues 
to the early formation process. These are 

- the outer boundary of the planetary system at 30 AU and 
inner boundary of the Kuiper belt at ~ 45 AU. 

- the much larger eccentricities and inclination angles of the 
Kuiper Belt objects, as a group, than the planetary orbits 

- the presence of precursors to the short-lived radioactive iso- 
topes ^''Fe and ^^Al in meteorites 

These properties strongly indicate that the birthplace of the 
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sun is likely to have been a cluster environment (Gaidos 1995 



IFernandez) |1997[ |Adams| |2010| . Perhaps the strongest con 
straint comes from the presence of the short-lived radioactive 
isotopes ^'^Fe and ^*A1 inferred from their precursor products 
in meteorites ( [Wasserburg et al.| |2006[ |Thrane et al.\ |2006[ 
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have half-lives of 2.62 Myr and 0.717Myr, respectively. Since 
they decay so quickly, their presence in meteorites means that 
they must have been incorporated into the early Solar System 
very soon after they were produced. 

A supernova is the most likely source for these short-lived 
radioactive isotopes, in particular the isotope ''"Fe is extremely 
difficult to produce otherwise than by stellar nucleosynthesis 
(Looney et al. 2006, Williams & Gaidos 2007, Gounelle et al. 
2009). A range of supernovae progenitor masses Msp are possi- 
ble, and no mass scale produces perfect abundances. However, 
these studies suggest that stars with M^, « 25M0 provide the 
best fit to the ensemble of short-lived radioactive nuclei (Adams 
2010, and references therein). Stars of such high mass usually 
form only in relatively massive clusters. External enrichment of 
short-lived radio-isotopes suggests a cluster with at least > 
4000 in order to have a reasonable chance of producing a star 
with Mjp=25Mo(Lee et al. 2008). 

The relatively high concentration of these short-lived nuclei 
means not only that they should have been formed just after 
the formation of the Sun but also fairly close to it. If not, the 
cross-section of the solar protoplanetary disc would have been 
too small to account for the amount of material. Several authors 
have estimated the maximum distance of the solar system when 
this supernova exploded, obtaining values between 0.2 and 2 pc 
(for a discussion, see Adams 2010). 

The outer boundary of our planetary system (30AU) and in- 
ner boundary of the Kuiper belt at 45 AU are likely significant 
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Table 1. Constraints on the properties of the solar birth cluster, 
where is the number of cluster members and pc is the central 
cluster density. 



variable 


value 


limiting factor 


N 


>4000 


chemical composition 


N 


< several 10"* 


radiation field ^ 


Pc 


> 10^ Mopc-3 


Sedna orbit ' 


Pc 


< lO'' Mopc--'' 


Sedna orbit * 



Note: 'Lee et al. (2008), ^ Adams (2010), ^Brasser et al. (2006), 
^Schwamb et al. (2010). 



for the early history of the solar system, too. This outer edge 
is relatively sharp - as best illustrated by the low mass con- 
tained in the total of the large number of Kuiper belt objects, 
which is only ~ 0.01 - 0.1 times the mass of the Earth (Bernstein 
et al. 2004). The protoplanetary disc from which the solar sys- 
tem eventually formed was most likely considerably larger at the 
time of planet formation. Possible processes that could have led 
to the truncation of the disc include; gravitational interactions 
with other cluster member s pda et al. |2000[ Kobayashi & Ida 



200T]|Pfalzner et al.j [2005] [Kenyon & Bromley 1, 12004', 'S purzem 
et al.| |2009)l, photcTevaporation by nearby massive stars (Adams 
et al. 2006"; Owen et al. 2010; Mitchell & Stewart, 20TTfor 
the supernova explosion itself (Chevalier 2000). Each of these 
scenarios require a cluster environment of relatively high stellar 
density. 

An alternative explanation for the sharp edge in the mass dis- 
tribution in the solar system would be the Nice model (Gomes 
et al. 2005), which suggests that the migration of the giant plan- 
ets caused a planetesimal clearing event leading to a late heavy 
bombardment (LHB) at 880 Myr. However, Booth et al. (2009) 
find that collisional processes are important in the Solar system 
before the LHB and that parameters for weak Kuiper belt objects 
are inconsistent with the Nice model interpretation of the LHB. 

The much larger eccentricities and inclination angles of the 
Kuiper Belt objects as a group compared to the planetary or- 
bits, also point to an encounter being the most likely cause. The 
most prominent example is Sedna with a periastron of 76AU, an 
aperastron of 937AU and an eccentricity of 0.8527. It is well be- 
yond the reach of the gas giants and could not be scattered into 
this highly eccentric orbit from interactions with Neptune alone 
( Gomes et al. 2005j). The most straightforward explanation for 
these high eccentricities would be some type of dynamical inter- 
actions (Morbidelli & Levison 2004) - either a close single en- 
counter or a wide binary solar companion (Matese et al. 2005). 

If an encounter is responsible for the disc cut-off and the 
Sedna orbit, then it must have taken place when the stellar den- 
sity was still relatively high. It is generally assumed that clusters 
expand significantly after gas expulsion which is accompanied 
by a dramatic drop in stellar densities. The time scale of this 
cluster expansion was until recently not well known. The en- 
counter must have occured early on in the history of the solar 
system (Malhotra 2008) as otherwise the dynamical interactions 
due to passing stars would have not only lead to the truncation of 
the disc and the high eccentrcities in the Kuiper Belt objects, but 
would also have perturbed the orbits of the planets, stripped the 
comets away and led to portions of the Inner Oort Cloud becom- 
ing unstable ( Gaidos , 1995] Adams and Laughhn |200 1). The 
fact that all of the planetary orbits are nearly in the same plane 
(inclination angles < 3.5°), their low orbital eccentricities(<0.2) 



and the presence of the Inner Oort Cloud all point to a lack of 
severe disruption after the solar system was fully formed. This 
provided an upper limit on the stellar density in the solar birth 
cluster for ages > 30Myr. 

Recently a number of authors tried to translate constraints 
derived from the meteorite composition, the disc cut-off and the 
Sedna orbit into properties of the solar birth environment. These 
are basically limitations on the number of stars in the solar 
birth cluster and its central density pc- For a summary of their 
results see Table 1 . These conditions do not necessarily need to 
be fulfilled over a long period of time. So the condition on 
based on the chemical composition needs only to be fullfilled 
before the supernova explodes, whereas the upper limit based on 
the radiation field needs to be fulfilled after gas expulsion from 
the cluster - before gas expulsion photo evaporation is much less 
efficient - until about 5 Myr, when discs are anyway largely dis- 
persed (Hernandez et al. 2007). The conditions on the cluster 
density must only be fulfilled at the instant of the encounter 

Most investigations conclude that a membership of ^ 10^ 
to several 10'* is the most likely scenario for the solar birth clus- 
ter. The lower limit of this estimate is based on the demand 
that the solar birth cluster should have contained enough stars 
in order to that there is a reasonable likelihood that it hosted 
the required supernova progenitor (for details see Adams 2010). 
Chemical considerations also suggest that the Sun formed in 
the presence of strong FUV radiation fields, where rough esti- 
mates indicate a birth cluster with A^ > 4000 (Lee et al. 2008). 
Really massive clusters with A^ >10^ have been considered as 
birth place of the solar system (Williams & Gaidos 2007, Hester 
et al. 2004) too. However, the radiation fields provided by such 
massive young clusters are potentially disruptive. Nevertheless, 
the early solar nebula could survive in clusters with A^ ~ sev- 
eral 10"* provided it spends enough time in the outer cluster re- 
gions(Scally & Clarke 2001, Mann & Williams 2009). 

The cut-off radius of the disc is also used to constrain the 
membership of the solar birth cluster The basic method here 
is to determine the required periastron distance rp^^i of the Sun 
with respect to the star encountered to obtain a cut-off in the 
disc. In a second step this is related to a certain cluster mem- 
bership. The lower li mit of rpgri is deri ved from the results of 



Hall et al.| ( |l996| l and [Kobayashi & Ida| ( [2001[ ), which state that 



parabolic, prograde encounters with solar-type stars lead to a re- 
duction of the disc size to approximately 1/2 - 1/3 of the perias- 
tron distance. Assuming a typical disc size of lOOAU, the cut-off 
radius translates into a periastron distance of 



pen 



lOOAU to 

a grazing encounter. Here a number of simplifying assumptions 
(parabolic, solar-type star) are made which might require recon- 
sidering (Olczak e t al] [20081 |20T 0). Other authors (Morbidelli[ 
and Levison 2004{ [Looney et al?) 200 6 ) come to somewhat dif- 
ferent periastron values. So to cover the entire parameter range, 
we consider periastra from 100 AU to lOOOAU to be potentially 
solar-system forming encounter events. This translates into a so- 
lar birth cluster in the approximate range A^ - IQr'- lO'* (Adams 
2010). 

To recap, the abundances of short-lived radioactive isotopes 
and the cut-off at 30AU both strongly hint at a solar birth cluster 
size in the range of 4000 < N < IQr' . Assuming that the stellar 
masses in the cluster follow the initial mass function (IMF) with 
an average stellar mass of niav - 0.5Mq this is equivalent to a 
total cluster mass rrid of 2000 Mo< ntd < 0.5 x IO^Mq. 

What about the stellar densities in these clusters? The den- 
sity of the solar system birth cluster is so far mainly constrained 
by considering what kind of encounters are necessary to result in 
a Sedna-like orbit. It is a point of debate whether there is an en- 
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counter scenario that can account for both features - cut-off and 
the Sedna orbit - resulting from a single encounter, or whether 
two separate encounters are necessary. 

Performing investigations of the origin of Sedna's orbit, 
Brasser et al. (2006) req uired central cluster d ensities higher 
than 10^ Mopc"^ whereas 



Schwamb et al. 



(2010 1 ruled out cen- 



tral densities >10 Mopc . Both worked with somewhat dif- 
ferent periastri than Adams(2010), which can be summarized in 
that an encounter between 100 - 1000 AU was required to result 
in the high eccentricities of the Kuiper belt objects. Brasser et al. 
(2006) and Schwamb et al. (2010) give limits to the central, and 
not the mean density of the cluster, because due to mass segrega- 
tion the most massive stars of a cluster are usually located close 
to the cluster center Since the solar system must have formed 
close to a supernova which had a massive star as progenitor, the 
Sun most likely resided close to the cluster center 

In summary, it can be said that the presence of short-lived 
isotopes, the fall-off of the mass distribution of the solar system 
and orbit of Sedna all combine to constrain the solar birth cluster 
to a membership of a few thousands to ten-thousands of stars 
with a central density of 10^ < Pc ^ 10^ Mopc"^. 

Recently it was found that massive clusters evolve in specific 
ways after gas expulsion (Pfalzner 2009). In this paper we inves- 
tigate how this knowledge can be applied to deduce the birth 
environment of the Solar system. Section 2 recaps the tempo- 
ral development of massive clusters. In Section 3 we show by 
means of simple approximations that the Solar system is un- 
likely to have developed in a star burst cluster. Section 4 shows 
the results of detailed modelling of the temporal developement 
of solar type star-disc systems in leaky clusters. The limitations 
and consequences of these results are discussed in Section 5. 



2. Cluster environments 

The requirement of the mass of the solar system birth cluster 
to be in the range 2000 Mq< nid < 5 x IO^'Mq puts it in the 
class of massive clusters. In the Milky Way massive young star 
clusters (t^. < 20 Myr, > 10"^ Mg) are observed to have mean 
densities ranging from less than 0.01 to several 10^ Mopc"^ and 
radii from O.lpc to several tens of pc. This means clusters with 
properties, in terms of mass and density, demanded for the solar 
birth cluster exist also at the present time in the Milky Way. 

Recently it was demonstrated that such massive clusters 
evolve in time only along one of two well-defined density-radius 
tracks - strongly suggesting a bi-modal cluster formation (see 
Fig.l) (Pfalzner 2009). The two cluster groups differ consider- 
ably in their properties: One group, termed starburst clusters, 
comprises young compact clusters which develop from a radius 
Rc of ~ 0. Ipc at 1 Myr to about l-3pc at 20 Myr. Simultaneously 
the density drops from the initially very high mean stellar density 
of >100^ Mopc"^ to 10 Mopc"^ during that timespan. All these 
starburst clusters are born with the same mean cluster mass, ra- 
dius and density and then simply diffuse as R^^ while more or 
less retaining their mass. The expansion basically proceeds lin- 
early in time, i.e. as R^. oc t^ with an expansion velocity of 0.1-0.2 
pc/Myr. 

The other group, termed leaky clusters, contains many clus- 
ters that are usually termed OB associations. This group consists 
of clusters covering the same age span, but having on average 
larger radii and lower mean cluster densities (<I - 10^ Mopc"^). 
Starting with cluster radii of ~2pc at 1 Myr, they also expand 
in a time sequence, however, here the cluster density decreases 
as R^"^ rather than R^-*. In these clusters in addition to diffusion, 
the clusters lose mass. The cluster size does not increase linearly 




Cluster radius[pc] 



Fig. 1. Observed cluster mass (filled symbols) and derived den- 
sity (open symbols) as a function of cluster size for clusters more 
massive than 10^* Mgas published by Pfalzner (201 1). The values 
were taken from Figer (2008), Wolff (2007) and Borissova et al. 
(2008) and references therein. 



with time, but approximately as t"-^ "'^ and reaches about 25pc 
at 20 Myr. 

This restriction of the development of clusters more mas- 
sive than IO^'Mq along only two possible developmental tracks 
has profound consequences on the formation history of our solar 
sytem: The solar system inevitably developed either in a star- 
burst or a leaky cluster The now available knowledge of how 
the average cluster density pc changes with cluster age tc in these 
two environments, reduces the scenarios of the solar birth clus- 
ter to only two possible developmental paths. It now remains to 
determine which of the two environments - starburst cluster or 
leaky cluster - is the more likely origin of the solar system. 

3. Solar system in a starburst cluster? 

Starburst clusters retain more or less their complete mass over 
the first 20 Myr. They are likely to remain a bound entity until the 
galactic tidal field disrupts them. This probably only happens for 
ages tsb >100Myr. A statistical argument against starburst clus- 
ters as solar birth hosts would be that probably only ~ 10% of the 
stellar population is born within systems that remain gravitation- 
ally bound over timescales longer than 100 Myr (e.g., Battinelli 
& Capuzzo-Dolcetta 1991, Adams & Myers 2001). Most stars 
are born within clusters that dissolve quickly, after only a few 
tens of Myr (Lada & Lada 2003, Porras et al. 2003). 

However, the strongest argument against starburst clusters 
is their very high initial density of p'^^ >10'' Mopc"^, combined 
with the high number of massive stars. From the density-size and 
size-age relations derived in Pfalzner (2009), it follows that the 
development of the mean stellar density psb in a starburst cluster 
environment can be approximated by (see also dashed line in 
Fig. 2) 



p,b ^ 5 X 105f;3±°-3[Mopc-3]. 



(1) 



where and in the following f^. is in units of Myr. Before com- 
paring this density development in starburst clusters with the re- 
quirements derived by Brasser et al. (2006) and Schwamb et al. 
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Fig. 2. Mean cluster density as function of the cluster age for 
starburst (circles - observed values and dashed line - interpola- 
tion as in Pfalzner 2009) and leaky clusters(triangles and solid 
line, observations and interpolation, respectively). The red box 
shows the density requirement for solar-system forming encoun- 
ters to be likely. The open squares show the parameters of the 
clusters modelled in chapter 4. The thin dashed line shows a fit 
through the data case where the Arches cluster is not considered 
in the starburst cluster sequence. Here only the errors in cluster 
age are given, the errors in cluster densities can be found in Fig. 
2 of Pfalzner 2009. 



(2010) for the solar birth cluster, one has to correct for the fact 
that the values by Brasser et al. and Schwamb et al. are cen- 
tral densities whereas Eq. 1 describes the mean density. The rea- 
son why they consider central cluster densities is that many star 
clusters exhibit mass segregation at early ages (< IMyr) with 
the most massive stars preferrentially being located close to the 
cluster center. From that it follows that most likely the Sun was 
located near the cluster center as otherwise it could not have been 
in close vicinity to the exploding supernova. 

The mass density values quoted in Pfalzner (2009) are ob- 
tained by dividing the cluster mass by its volume - implying an 
average density. The stellar density in a cluster is in fact a strong 
function of the radial distance to the cluster center. Consequently 
the central density of a cluster is usually much higher than its 
mean density. For model clusters of the Orion Nebula Cluster 
(details of the simulation method can be found in Olczak et al. 
2010) we find - depending on the definition of 'central' - the 
central cluster density ranges from a factor 5-100 times higher 
than the average stellar density obtained at the half-mass radius. 
Observations confirm this fact at least for low-mass clusters: for 
example, Guthermuth et al.(2003) found that in all of the three 
embedded clusters that they observed, the mean volume densi- 
ties were of the order 10^-10^* Mopc^^*, whereas their peak vol- 
ume densities ranged from lO'^-lO^ Mopc^^*. Whether this holds 
as well for high-mass clusters needs further investigation. 

So the requirement by Brasser et al. (2006) and Schwamb 
et al. (2010) that the central density in the solar system birth 
cluster should He in the range lO^Mopc"-' < p** <10^Mopc"-' 
translates into a mean density of lOMopc""* < p;„ < lO'^Mopc"^. 
In Fig. 2 this area is depicted by a red box. It can be seen that the 
density requirement for the solar birth cluster overlaps with the 
development of starburst clusters at ages > 4 Myr 

However, at younger ages starburst clusters have an up to 10^ 
times higher density. At such high densities many of the stars 
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Fig. 3. Probability for encounters closer than lOOAU taking 
place as function of cluster age. The solid and the dashed lines 
show the estimates for starburst and leaky clusters, respectively. 
The thick and thin lines distinguish between the entire cluster 
and the central area only. 

would be completely stripped of their discs due to gravitational 
interactions ( |01czak et al.||2010l l. 

The average time Te„c until a star undergoes an encounter 
with a certain periastron distance rp^ri can be approximated by 
( |Davies[[20n] ) 



33Myr 



lOOpc 



n 



Ikm/s 



10^AU\/M, 



10 



(2) 



where n is the stellar number density in pc 



a the mean 

relative velocity at infinity of the cluster stars and m the mass 
of the considered star in Mq. With the cluster mass x Q.5Nc 
Mq, where Nc is the number of cluster members, it follows that 
P[MqPc^^] « 0.5n. Therefore the encounter probability in a star 
burst cluster can be approximated by 



'200 Mp 



T7nc ^ 33Myr|— — j-p,i, 



Ikm/s 



10'AU\/M, 



' pen 



_i0 

m 



(3) 



In a simplified picture an encounter closer than lOOAU 
would lead to a truncation of the disc to r^i^c <30AU (Adams 
2010). This deprives the disc of the material necessary to form 
the outer planets, so the resulting planetary system would no 
longer resemble our solar system. Equation 3 can be used to es- 
timate the likelihood of such a destructive encounter to occur in 
a starburst environment. Fig. 3 (solid lines) shows the likelihood 
of a destructive encounter (r^^,., <100AU) in a starburst cluster as 
a function of cluster age using the density development given by 
Eq. 1 . The thin line shows the probability for such an encounter 
to happen in the entire cluster and the thick line the result for the 
cluster center, assuming the density there to be on average about 
20 times denser then the mean cluster density. Here disc destruc- 
tion is assumed if each star is likely to experiences at least one 
encounter closer than 100 AU during a 1 Myr time interval. Or in 
other words, the destruction probability is 1, if t**^. < IMyr and 
l/T**^[Myr] for T^J^^ > IMyr This definition does not include re- 
peated encounters. So in reality disc destruction would be even 
higher than indicated in fig. 3. 
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In both cases the cluster environment is extremely collisional 
and encounters closer than lOOAU are very common. For the 
outskirts of the cluster this changes after 3 Myr but for the central 
area continues up to 7-8Myr. Disc destruction happens on fairly 
short timescales. For example, for a solar-type star in a cluster 
with pc >10^Mq pc"^, which is typical for starburst cluster envi- 
ronments younger than < 2Myr, an encounter closer than lOOAU 
can be expected already after t* <s;0.3Myr. The disc destruction 
happens on such short timescales (<lMyr) that it would be dif- 
ficult to form a planetary system before the disc is destroyed. 
Though one can not completely exclude the possibility of so- 
lar system having developed in a starburst cluster environment, 
it seems rather unlikely. Moreover, photo-evaporation ( Adarns] 
|2010| l, which has not be included in this estimate, would lead to 
additional disc destruction. 

Recent observations by Stolte et al. (2010) suggest that there 
exist discs around some of the stars in the starburst cluster 
Arches, however, the vast majority is devoid of inner discs. They 
only find discs around B-type stars with a frequency of 6% + 2%. 
For solar-type stars the likelihood of retaining a disc would be 
even lower (Pfalzner et al. 2006). Assuming for a moment that 
the solar system happened to develop in one of such rare cases 
of its disc surviving the first Myr, the chances of discs surviving 
even longer would actually increase, because starburst clusters 
become much less dense as they age. Nonetheless, at an age of 
20 Myrs the average stellar density would still be several 100 
Mopc"-'. This is such a high density that, even if one assumes 
that some discs in a starburst environment survive long enough to 
form a planetary system, it would inevitably lead to a strong per- 
turbation when the planetary system is close to being completed. 
Interactions with other cluster members would inevitably lead to 
disturbances destroying the near circularity and co-planarity of 
the solar systems. 

The only way that our solar system could have formed as part 
of a starburst cluster would be if it had formed at the outer edge 
of this cluster, always remained there and never transversed the 
cluster center. A rough estimate using Eq. 2 would result in the 
requirement that the solar system would then have to have been 
located at 4-5 times the half-mass radius at the time of cluster 
formation. 

It has been suggested that Arches might be exceptional for 
Milky Way clusters as it is thought to evolve quite differently 
from the clusters in the disk because of the unusually strong 
tidal fields (only ~30 pc away from the Galactic center). If we 
exclude this cluster from our investigation, this means that the 
size of the cluster at the time of gas expulsion could be larger 
than anticipated by including it. The thin dashed line in Fig. 2 
shows the density developemnt excluding Arches. The cluster 
expansion starts at lower densities on the developmental track in 
Fig. 1 . However, even in this case the cluster would go through a 
phase of 3-4 Myr where basically every star in the cluster center 
undergoes at least one encounter closer than 100 AU encounter 
in every 1 Myr timespan. So it still holds that nearly all discs 
are likely destroyed in the cluster center However, there is an 
increased change that discs could survive in the cluster outskirts. 
The requirement reduces to the condition that the solar system 
would then have to have been located at 2-4 times the half-mass 
radius at the time of cluster formation. 

However, at the same time a massive star would have to 
have been at these outskirts of the cluster to account for the ob- 
served chemical composition. Although this scenario cannot be 
excluded, it seems fairly unlikely. In the following we will see 
that it is much more likely that the solar system developed in 
leaky cluster environment. 



4. Solar system in a leaky cluster 

For the leaky clusters the age-dependence of the average cluster 
density p/ can be approximated by (Pfalzner 2009 | l 



Pi 



150f;2''*°-2[MQpc-3] 



for t> 1 Myr 



(4) 



where is in units of Myr. This is illustrated by the solid line in 
Fig. 2. Early on in the cluster development the mean cluster den- 
sity will be approximately a few times 10 Mo/pc-' to lOOMo/pc^ 
and will rapidly decrease over the following 20 Myr. Looking 
again at Fig. 3 (dashed lines) one sees that disc destruction ap- 
proximated from 

is not an issue in the leaky cluster environment. Even early on in 
the cluster development, less than 10% of the stars in the cluster 
center are likely to undergo encounters closer than lOOAU and 
in the entire cluster it < 1 % of solar-type stars that are afifected. 

Fig. 2 shows that the density requirement for the solar birth 
cluster overlaps with the development of leaky clusters for ages 
tc < 4 Myr At later times the cluster density is considerably 
lower making a solar system forming encounter unlikely. This is 
the central result of this study: with a very high probability, the 
solar system developed in a leaky cluster. 

This development of the stellar density in such leaky clus- 
ters also neatly solves the seemingly contradictory situation that 
a dense birth cluster is required to provide for the nuclear en- 
richment and explain Sedna's orbit, but at the same time, a less 
interactive environment is needed to avoid disruptive dynamical 
interactions with the formed solar system. 

In the following we want to go beyond the simple estimates 
of Eqs. 2-4 and determine what the development of the solar sys- 
tem in a leaky cluster means in terms of its history by performing 
numerical simulations of its encounter dynamics. 



4.1. Numerical method and initial conditions 

The clusters identified in Fig.l as part of the leaky cluster se- 
quence are all exposed clusters, where the gas has been expelled 
from the cluster Although gas expulsion seems to be a vital in- 
gredient for the cluster expanding in the observed way, there are 
currently many unknowns setting the course of the following 
cluster dynamics and mass loss in the leaky cluster sequence. 
These uncertainites concern the time scales of the gas expulsion, 
the dynamical state - equilibrium, sub -or supervirial - before gas 
expulsion and to what degree early substructuring plays a role. 

Therefore instead of trying to include these unknown pro- 
cesses here, existing simulations (Olczak, Pfalzner & Eckart 
2010) of clusters in virial equilibrium at different densities are 
used to obtain a first estimate of the dynamics of the early so- 
lar environment as a function of cluster age. The densities of the 
modelled clusters are indicated in Fig. 2 as open squares. The 
densities of these model clusters decrease very slightly over the 
considered timespan of 1 Myr This means the observed cluster 
density development is replaced by snapshots of clusters of fairly 
constant density. The actual numerical method and the model set 
up is described in detail in Olczak, Pfalzner & Eckart (2010). 

The main aspects of the model can be summarized in the 
following way: All simulations were performed using the direct 
many-body code Nbody6-n- (Aarseth 2003; Spurzem 1999). For 
simplicity it was assumed that all stars are initially single and no 
primordial binaries are considered. This seems to be justified. 
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as Adams et al. (2006) and Brasser et al. (2006) found in their 
investigation of potential solar system environments that binary 
scattering events do not influence the results significantly. The 
effect of stellar evolution has not been included. This is justified 
by us modelling only the first 5 Myr of the cluster development 
in detail, where stellar evolution plays a minor role. 

The stellar masses in the cluster were sampled from the ini- 
tial mass function given by Kroupa(2001). The cluster members 
cover the mass range from the lower mass end of O.OSM© to 
the currently accepted upper mass limit of 150 M0(Koen 2006; 
Oey & Clai-ke 2005; Mai'z Apellaniz et al. 2007; Zinnecker & 
Yorke 2007; Weidner & Ki'oupa 2006; Figer 2005). The velocity 
distribution is Maxwellian and the cluster is assumed to be in 
virial equilibrium. In reality leaky clusters at that phase of their 
development are actually supervirial as the cluster still adjusts 
after the gas explusion process. Bastian & Goodwin (2006) find 
for somewhat larger clusters than we study here and with star 
formation efficiencies of x 30%, that it takes several 10^ years 
until the clusters reach their new equilibrium state. However, not 
really knowing the gas expulsion dynamics and star formation 
efficencies in leaky clusters in detail, we resort to the easiest as- 
sumption namely treating the virial equilibrium case. We will 
discuss later in how far this might influence our results. 

The cluster density profile was chosen in such a way that it 
resembles the observed ONC profile (McCaughrean et al.. 2002| l 
after IMyr of simulation. We chose an ONC-like density pro- 
file because it is the best known massive cluster just before gas 
expulsion. As we model the early stages of exposed clusters we 
assume that the shape of the density profile does not alter sig- 
nificantly (at least in the central areas) over the timespan of the 
first a! 5 Myr. From much older open and globular clusters it 
is known that they have less centrally condensed Plummer-type 
profiles, but this development happens on timesclaes of >100 
Myr to Gyr. 

Due to the adopted stellar number density distribution, which 
is roughly represented by p = por^' , the density of the mod- 

els scales as the stellar number, N - p(r)r~drdD. cc poR.. 
Different initial stellar densities are obtained by varying the stel- 
lar numbers to 1000, 2000, 4000 (ONC), 8000, 16 000, and 32 
000 (see Table 2). These density-scaled cluster models have been 
simulated with the same initial size of R=2.5 pc corresponding 
to the size of the Orion Nebula Cluster (ONC). The central den- 
sity in the different models covers the range of a!l-42xl0'*pc"^. 

For each cluster model a set of simulations has been per- 
formed with varying random configurations of positions, veloc- 
ities, and masses, according to the given distributions, to min- 
imise the effect of statistical uncertainties. For the clusters with 
1000, 2000, 4000, 8000, 16 000, and 32 000 particles, a number 
of 200, 100, 100, 50, 20, and 20 simulations, respectively, were 
performed to ensure statistical robustness. Due to the relatively 
high particle number in the simulation the encounter statistics 
of the stars close the center is not influenced by the change in 
simulation particles (see Olczak et al. 2010). Therefore we can 
treat the obtained encounter probabilities as direct function of 
the cluster densities. 

By contrast to previous studies we do not consider the simu- 
lation of clusters of different densities as different possible solar 
system birth environments, but as consecutive stages in the de- 
velopment of a leaky cluster environment which we can assign 
to certain cluster ages. In the following this is described in more 
detail. 



Table 2. Cluster models, resulting mean cluster densities and 
age, where these densities are reached in the leaky cluster de- 
velopment. The last column shows the probability for a solar- 
type star located within 0.7 pc from the center to ungergo an 
encounter between lOOAU and lOOOAU 



No. of stars 


mean density 


cluster age 


encounter 




Mo/pc3 


Myr 


probability 


32000 


4.2 X 10^ 






16000 


2.1 X 10"* 




634% 


8000 


1.0 X 10* 




41.9% 


4000 


5.3 X 10' 


1.8 Myr 


28,4% 


2000 


2.7 X 10' 


2.3 Myr 


21,3% 


1000 


1.3 X 10' 


2.8 Myr 


16,9% 



100% 



80% 



o 

u 

3 
a 



60% 



40% 



20'^ 



B- □ within 0.7 pc 
• — • within 0.5 pc 
* ■ * < lOOAU 




10 15 20 

3 -3 

mean Stellar density (10 pc ) 



25 



Fig. 4. Percentage of solar-type stars (0.8Mq< Mj < 1.2Mq) lo- 
cated within 0.7pc (solid line) or 0.5pc (dashed line) from the 
cluster center at the moment of the encounter which has a peri- 
astron in the range lOOAU -lOOOAU during the first Myr of the 
cluster development as a function of mean cluster density. The 
mean density is the spatial and temporale mean during this 1 
Myr timespan. The dotted line shows the percentage of stars that 
have an encounter closer than 100 AU. 



4.2. Encounter statistics 

While simulating the stellar dynamics of above described clus- 
ters, for solar-type stars the parameters (mass ratio, periastron, 
eccentricity) of every encounter closer 1000 AU is recorded. 
Here we define solar-type stars as a stars in the mass-range 0.8- 
1.2 Mo . 

The Sun must have been close to a massive star that exploded 
as a supernova and enriched the solar system with ^''Fe and ^^Al. 
To consider this we again take the ONC as a guide for the typ- 
ical cluster environment. Here the result of mass segregation is 
that most massive stars are located within 0.5 pc of the cluster 
center. To model the early history of the Solar system we there- 
fore consider only stars within 0.5pc and 0.7pc form the cluster 
center. The latter denotes the central cluster area of 0.5 pc plus 
the estimated distance of the Sun to the supernova of 0.2pc. 

Figure |4] shows the probability of a solar-type star located 
within 0.5pc (solid line) and 0.7pc (dashed line) from the cluster 
center having an encounter between 100 AU and lOOOAU dur- 
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Fig. 5. The dots show the simulation results of the percentage 
of solar-type stars (0.8Mq< M, < 1.2Mq) located within 0.5pc 
from the cluster center which have a single encounter with a pe- 
riastron in the range lOOAU -lOOOAU as a function of cluster age 
for leaky cluster environments. The triangles shows the percent- 
age of encounters closer than 100 AU. The solid and the dashed 
lines show the equivalent estimates according to Eq.2. 



ing 1 Myr of cluster developement as a function of the cluster 
density. Stars that have encounters closer than 100 AU before or 
after in addition are excluded, because such encounters would 
truncate the disc to such a degree that they would be too small 
for developing into a solar-type planetary system. As expected 
the likelihood of encounters between 100 and lOOOAU increases 
considerably with higher cluster density. However, so does the 
value of encounters closer than lOOAU - the proportion of "de- 
stroyed" discs (dotted line in Fig 4] ). Therefore in very dense 
clusters (p > 2x lO'*) the likelihood of a solar system forming en- 
counter decreases again, because an increasing number of discs 
have an encounter closer than lOOAU that makes it unsuitable 
for developing into a solar system. 

Given the encounter likelihood at different central densities, 
we can now deduce the encounter statistcs during the develop- 
ment in a leaky cluster with age using eq. 3. Table 2 shows how 
the different densities correspond to ages in the leaky cluster se- 
quence. No densities above ~ 5x 10^ Mopc^^* are reached in 
the leaky cluster sequence. However, this does not necessarily 
mean that the cluster never went through a phase of higher stel- 
lar density. Such a situation might have existed in the earlier 
embedded phase while the cluster was still in its formation pro- 
cess. However, little information about the stellar densities in 
forming massive clusters exist and existing knowledge can be 
interpreted very differently (for different views see for example 
Pfalzner 2011 and Parmentier et al. (subm.)). Therefore we re- 
strain this study to the exposed phase only. 

Figure 5 shows the probability for an encounter at a distance 
of 100 AU to 1000 AU in such an evolving leaky cluster. Here 
table 2 and eq. 4 have been used to equate a given density to a 
cluster age. The encounter rate is obtained as an average over 
0.5 Myr During that time interval the cluster density decreases 
faster in the leaky cluster development than in the simulated sys- 
tems, which are in virial equilibrium. Therefore the encounter 
rates have to be regarded as upper limits which is reflected by 
the larger error bars towards lower encounter rates. 



Fig. 6. Distribution of the perturber mass of the encounters be- 
tween 100 AU and 1000 AU for solar-type stars in leaky clusters 
at 1.8 (dashed line) and 2.8 Myr (solid line). 



Figure 5 shows that in a leaky cluster environment solar-type 
stars have a probability of ~ 30% of experiencing a potentially 
solar system forming encounter early on in the development. As 
the cluster density decreases with cluster age, so does the en- 
counter probability. This plot demonstrates that the highest like- 
lihood for a solar system forming encounter is early on in the 
cluster development (<2Myr). After 5 Myr of cluster develop- 
ment encounters below 1000 AU become extremely rare because 
the density in the cluster has decreased so much. 

It can be seen that the encounter rate measured in the N-body 
simulations (circles) is higher than the one expected from the 
simple estimate of eq. 5 (solid line). This is because there is grav- 
itational focussing around the massive stars, which themselves 
are concentrated in the cluster center (jPfalzner et al. 2006 Fujii 
|& Portegies Zwart 201 1), which in turn leads to a significantly 
higher number of encounters. The combination of gravitational 
focussing and mass segregation is not considered in the estimate 
of eq. 5. 

These more detailed studies support the idea that disc de- 
struction is less of a problem in leaky clusters. Although the 
simulation value of disc destruction (triangles) is higher than the 
estimated one (dashed line), still very few solar-type stars have 
an encounter that is closer than the threshold of 100 AU. 

Figure 6 shows the mass distribution of the encounter part- 
ners in a timespan of 1 Myr in these encounters between lOOAU 
and 1000 AU at an cluster age of 1.8 Myr and 2.8 Myr, respec- 
tively. It can be seen that there are two maxima in this distribu- 
tion - one for low-mass and one for high-mass perturbers. Most 
encounters are with low-mass stars ( 1112 < 0.8Mq), this is sim- 
ply due to the fact that they are the most common in the cluster. 
However, a considerable amount of encounters take as well place 
with the highest mass stars (> 10 Mg) of the cluster As Fig. 6 
shows it is actually just as likely for the Sun to have had an en- 
counter between 100 AU and lOOOAU with a solar-type star as 
with a high-mass star despite low numberof the latter. The rea- 
son is gravitational focussing. This is most pronounced later in 
the cluster development when the cluster density is lower This 
result shows that the often applied assumption that the perturber 
was also a solar-type star has to be reconsidered in future work. 
These new results do not even exclude that the supernova pro- 
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Fig. 7. Eccentricity distribution of the encounters between 100 
AU and 1000 AU for solar-type stars in leaky clusters at 1.8 
(dashed line) and 2.8 Myr (solid line). 



genitor itself might have been the encounter partner that led to 
the shape of the solar system. 

There is another point about the solar system forming en- 
counter that might require revision: the eccentricity. Fig. 7 shows 
the distribution of the eccentricities of the encounters for two 
different cluster ages as the leaky cluster environment develops. 
In the early cluster development strongly hyperbolic encounters 
are more likely. Hyperbolic encounters lead to considerably less 
mass loss in the disc than the equivalent do parabolic ones be- 
cause the interaction time for disc mass removal is much shorter 
(see, for example, Pfalzner 2004, Pfalzner et al 2005, Olczak et 
al. 2010). Only at ages >2Myr parabolic encounters start to hap- 
pen more often. 

The fact that highly hyperbolic encounters dominate in dense 
cluster environments might increase the chance that the solar 
system formed in a starburst cluster environment because the 
disc destruction timescale in these environments might increase. 
First indications that this might be the case were found by 
Olczak et al. 2012. However, the problem of the strong radia- 
tion field in starburst clusters remains. 

In summary, it can be said that the assumption that the per- 
turber that formed the shape of the solar system was a solar-type 
star on a parabolic orbit is highly questionable (see as well Dukes 
& Krumholz 2011). Future investigations will have to take the 
entire mass spectrum and a much wider eccentricity range into 
account to narrow down the encounter sceanrio that lead to the 
formation of the solar system. 

5. Discussion and conclusions 

Previous work puts the cluster population limit of the birth 
environment of the Sun at 10^ to 10^ members. Recently it 
was shown that such massive clusters only develop along one 
of two existing tracks in the density-radius plane - as either 
starburst clusters or leaky clusters(OB associations). As a 
consequence the solar system must have developed in one of 
these two environments. Using the now known temporal density 
development of starburst clusters, it has been demonstrated here 
that in starburst clusters the high initial central stellar density 
(>10^ Mq/pc^) despite its rapid decline leads almost inevitably 



to disc truncation over the first 5-lOMyr. The resulting discs are 
often too small to account for the geometry of the solar system. 
In addition, the high concentration of massive - and therefore 
very luminous - stars drives additional disc destruction. The 
combination of both effects, makes it unlikely that the solar 
system formed in a starburst cluster. 

The central results are: 

- The solar system most likely formed in a leaky cluster envi- 
ronment. 

- As a consequence, the solar birth cluster dispersed to a large 
degree over timescales of ~ 20 Myrs, where the cluster ra- 
dius increased approximately like r cc f^^ and the density in 
the cluster diminished as pilM^pc^^] = ISOf^T^'^'^"'^ leaving 
behind a cluster that only has at most 10-20% of its initial 
mass. 

- Due to the rapid decrease in density in the solar birth clus- 
ter, encounters played a role in shaping the solar system 
only very early on in the cluster development. At 1.8 Myr 
the probability for an encounter at an distance of 100 to 
lOOOAU is ^ 30%. The encounter most likely took place 
when the cluster was less than 5Myr old, because by then 
the likelihood for such an encounter dropped well below 
10%.. If the development to less concentrates profiles hap- 
pens on timescales <10 Myr, than we would overestimate the 
encounter frequency in the central parts at the latter stages 
(>3Myr) in our simulations. So our results at >3Myr can be 
regarded as upper limits. This means that cluster older than 
5 Myr would be even less likely to experience a solar system 
destroying encounter. 

- Hyperbolic encounters with low mass stars and parabolic en- 
counters with high-mass stars are the most common events 
for solar-type stars in solar system birth clusters. The former 
mostly happen in the very early stages (<2Myr), the latter 
preferrentially in the timespan of 2-4Myr. 

Perturbation of the already formed planetary system 
([Maknberg eTaT] [20071 PTT] [Lestrade et al .| |20TT] [Parker & 
Qu anz|[2011| l by encounters did not take place, because having 
formed in a leaky cluster the encounter probability for the solar 
system rapidly dropped as the cluster density decreased precipi- 
tately with cluster age. Even if the solar system was part of the 
remnant cluster after gas expulsion the stellar density would be 
< 1 Mopc"^, far to low to make a close encounter likely. 

For the same reason the presence of so many massive stars 
(10-100) in the solar birth cluster did not lead to its destruc- 
tion. Although the most massive stars only need 3-5 Myr from 
their formation to their explosion, even if these supemovae were 
closely packed at the start of the cluster development due to mass 
segregation, the rapid cluster expansion of the cluster can ac- 
count for the solar system being spared further disruption by su- 
pemovae at later times. 

The here investigated encounters within the distance range 
of 100 to 1000 AU lead not only to a cut-off in the disc but as 
well to a steepening of the mass distribution inside the remaining 
disc. Steinhausen et al (2011) showed that this could result in a 
r^^ -dependence as required for the minimum mass solar nebula 
(Desch 2007). 

In the still embedded phase the stellar density in leaky clus- 
ters were still perhaps somewhat higher (Pfalzner 2011, for an 
alternative view see Parmentier et al (subm)). If this were the 
case, possibly the encounter frequency in the embedded phase 
could have been higher than in the exposed phase. However, star 
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formation proceeds most likely in an accelerated way with most 
stars forming just before gas expulsion. In this stage the stellar 
density is only slightly higher than at the start of the exposed 
phase, so the sun most likely spent only a short time in this 
slightly denser environment. 

In this study we took the earlier results of the most likely 
periastron of the solar system-forming encounter having been 
between 100 and 1000 AU at face value. These former stud- 
ies often considered only encounters between two solar-type 
stars and/or parabolic orbits. However, the results obtained here 
strongly suggest that these assumptions are not necessarily justi- 
fied. Future work should consider as well hyperbolic encounters 
and include the entire mass spectrum of encounter partners. This 
might be especially important in the evaluation of the likelihood 
of the solar system having formed in a starburst cluster 

In the current study we included the cluster density devel- 
opment in a simplified way, by looking at simulations of clus- 
ters of different densities in virial equilibrium. Obviously leaky 
clusters are at that stage in their development not in virial equi- 
librium, most prominently illustrated by their mass loss. So in 
future studies it will be essential to include the gas expulsion 
phase to treat the cluster expansion in an appropriate way. 

We have recently modelled the cluster expansion of starburst 
and leaky clusters after gas expulsion. Preliminary results seem 
to show that the observational data are surprisingly well repro- 
duced by the simulations (Pfalzner & Kacmarek, in preparation). 
However, the fact that bound and unbound stars with the latter 
having a preferrential direction occupying the same space might 
influence the encounter statistics. This needs to be investigated 
in future work. In addition, we find indications that the leaky 
cluster sequence might only be representative for the develop- 
ment of the most massive clusters in the Galaxy (> 10^ Mq). 
Lower mass clusters, although starting with approximatedly the 
same cluster size, develop along a slightly different path, quickly 
decreasing so much in density that they fall below the detection 
limit. If these results are confirmed, the cluster mass range 2-10 
X IO-'Mq would require to be investigated in the light of these 
new results again. However, given the fast decrease in density 
in these lower mass clusters the time span of a high encounter 
probability is extremely short - probably in the range of IMyr or 
even less. 
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